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A SKEW NORMAL MODEL OF THE DOSE-EFFECT RELATION 

IN PHARMACOLOGY 

FRANCIS LAVERGNE AND NICOLAS MARIE 


Abstract. This paper deals with a skew-normal model of the relation be¬ 
tween a dose d > 0 and a quantitative measure E(d) of an effect of the ad¬ 
ministered drug. Precisely, E(d) is a measure of the therapeutic response or a 
measure of a side-effect. Some existing and additional properties of the logis¬ 
tic functions are proved, and a skew-normal model of the escape time of rats 
under an experimental antidepressant medication is provided. 
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1. Introduction 

This paper deals with a skew-normal model of the relation between a dose d > 0 
and a quantitative measure E(<J) of an effect of the administered drug. Precisely, 
E (d) is a measure of the therapeutic response or a measure of a side-effect. 

As explained in Prentice [10] and Brown j4], it is usual to put 

E (d) := / f{x)dx 

Jo 

where, fi £ R, a > 0 and / : R. —> R. is a probability density function. Often, / 
is a probit or a logit function. The parameters /i and a are estimated from the 
observations. For an application on clinical datas, for instance, see Verlato et al. 

EH- 

The skew-normal distribution has been already used to model the quantal response. 
In [15], Section 3.5, Wagner assumes that / is a skew-normal density function. On 
the skew-normal distribution, see Azzalini [5] and Chen et al. [5]. The basics on 
the skew-normal distribution are stated at Section 3.1. On the multidimensional 
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skew-normal distributions and an application in neurotoxicology, see T. Baghfalaki 
et al. j3]. 


The model studied in this paper is a family of random variables E := {E(d), d £ R + } 
satisfying the following conditions : 


(1) 

( 2 ) 

( 3 ) 


For every d £ M + , E(d) Af(£(d), w 2 (d), a(d)) where £, a £ C°(M + ,R) 
and u> £ C°(R+,]0, cso[). 

[i : d £ R + i—>• E[E(d)] is a logistic function. 

There exists do € R+ such that cr : d £ M + i —> -y/var[E(d)] is decreasing 
on [do,c»[, and 

lim a(d ) = 0. 

d—> <50 


Assume that E(d) is a measure of the therapeutic response of the administered drug 
at the dose d £ R+. It is usual to assume that E(d) is a Gaussian random variable. 
The skewness of the empirical distributions are taken into account in the model E 
of the therapeutic response in order to refine the choice of an optimal dose. Indeed, 
the skewness of the distribution of E(d) indicates if the therapeutic responses of the 
major part of the patients are over (positive skewness coefficient) or under (nega¬ 
tive skewness coefficient) the mean therapeutic response. Therefore, if the mean of 
E(d) is high enough, its skewness coefficient is positive and its standard deviation 
is small enough, then the dose d is admissible. In fact, ideally, the optimal dose 
should maximise the mean and the skewness coefficient, and minimize the standard 
deviation of the therapeutic response. 


Section 2 deals with some existing and additional properties of the logistic func¬ 
tions. Section 2 provides an approximation method of the parameters of the logistic 
functions in the most general case, which is used at Section 4. The proofs of the 
results stated at Section 2 are detailed at Appendix A. 


Section 3 deals with the estimation of the functions y and a, and then of the 
functions £, ui and a. So, for each admissible dose d > 0, the model E allows to 
simulate the measure of the effect E(d) multiple times for a better evaluation at the 
dose d. As mentioned above, the basics on the skew-normal distribution are stated 
at Section 3.1. 


Section 4 deals with a model of the escape time of rats under an experimental an¬ 
tidepressant medication. In psychopharmacology, since the medications are often 
administered for several months or years, it is crucial to find the smallest efficient 
dose. See Lavergne and Jay [ 5 ] about the efficiency of antidepressant medications 
for small doses. 

In j7], Holford and Sheiner studied the relationship between the therapeutic re¬ 
sponse and the elimination process of the drug. In a forthcoming work, the dose- 
effect model studied in this paper will be related to the fractional pharmacokinetics 
model studied in Marie [5]. 

2. Approximation of the parameters of the logistic functions 

This section deals with some existing and additional properties of the logistic 
functions. 


Consider n £ N* observations yi ,..., y n > 0 of a logistic function / at x ±,..., x n 
respectively. Many authors approximate the parameters of the function / by linear 
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regression on (x\, Z\),..., (x n , z n ), where 



; Vie {l,...,n} 


and 


y* := max y t . 

£E{l,...,n} 


The major drawback of that method is to assume that 

lim f(x) = y* or lim f(x) = y*. 

x — y —oo x—y oo 


This section provides an approximation method of the parameters of the logistic 
functions in the most general case, which is used at Section 4. The proofs of the 
results stated in this section are detailed at Appendix A. 


Definition 2.1. / : 


/( x) := h + 


is a logistic function if and only if 
1 


(Z 2 - M -1 + e mx +P 
with m G R*, j)£R, and 1 1 ,12 G 1 such that I 2 > l\. 


; Vcc G 


Throughout this section, / is the logistic function defined at Definition 2.1 


Proposition 2.2. The logistic function f satisfies the following properties : 

(1) If in < 0 (resp. m > 0), then 

lim f(x) = l\ and lim f(x) = I 2 (resp. lim f{x) = I 2 and lim f{x) = l\). 


x — y —00 


x—>—oo 


(2) If m < 0 (resp. m > 0), then f is increasing (resp. decreasing) on R. 

(3) The graph of the function f has a unique inflection point, at 

0 , = log(Z 2 - h) +P 

m 

Corollary 2.3. The parameter li is a solution of the following equation : 

2 0f'(0) _ n 

m - h ' 

Moreover, 


(1) 


log 


J(0)~f( 0) 

" /(0)-ii + 


h = 2 f(S) -l,, 771= - jtL andp = log 
J{d) - ‘1 


m-h 


i 

2[/(0)-*iL ' 


Proposition 2.4. The logistic function f is the solution of the following ordinary 
differential equation : 


( 2 ) 


y{x) = /(0) - TO / (y(u) - h) 


1 - 


h — h 


[y(u) - h 


du. 


Consider n G N* observations y \,..., y n of the logistic function / at x\ = 
0, X 2 ■ ■ ■, x n respectively. The end of the section is devoted to some methods to 
approximate the parameters to, p, l\ and I 2 of the logistic function /. 


On one hand, assume that the values of the parameters l\ and Z 2 are known. Let 
$ :]l 1 , Z 2 [—>• R be the map defined by : 


<h(:r) := log 


1 


x — l\ 


1 


I 2 — k 


; Vx G]Zi, Z 2 [- 
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By Proposition |2.4| together with the change of variable formula for the Riemann- 
Stieljes integral : 


$[/(*)] 


aw)] + r <s>'[f(u)]df( U ) 
*[/(o)]- ° 


m[ &[f(u)][f{u)-l i] 
Jo 


[/(«)-«!] 

h ~ n 


du 


p — mi. 


So, if l\ and I 2 are known, the parameters m and p can be approximated by linear 
regression on (aq, Zi),..., (x n , z n ), where 

Zi ■= ®(Vi) i Vi G n}. 


Precisely, 


n — 


1 (E*J +E : 


m n := -- 


0=1 / \i=l 


and 


Pn := 


A n x-v. n 

1 V— win 

T / . ^ T / , x i 

n — 1 n — 1 

4=1 4=1 


are some unbiased estimators of m and p respectively. This method is classic (see 
R.F. Gunst and R.L. Mason [5]). 


On the other hand, assume that the values of the parameters l\ or I 2 are unknown. 
In that case, the transformation 4) cannot by applied to y \,..., y n . In the current 
paper, an alternative method of approximation is provided. 


Since card({l ,...,n — 1}) < 00 , the maximization problem 


max 

2 E{l,...,n— 1} 


Vi +1 - Vi 

%i -\-1 


has a unique solution n(9). So, 


0 n ■= 


x n(8)+l + x n(8) 


yn{6 )+1 + Vn(6) 
2 


and 8 n := 


2/n(0) +1 Vn(0) 
%n(6 )+1 3Cn{6) 


define some converging approximations of 0, f(0) and f'{6) respectively. 


If the value of l\ is known but not the value of I 2 , then by Corollary |2.3| : 

h{n) := 2y n - l x 

defines an approximation of Z 2 - So, the previous method is adaptable by replacing 
by the map <!>„ :]Zi, > K defined by : 


®n( x ) ■ = log 


1 


1 


X — l\ 


Vx €]li,h(n)[. 


l 2 (n) - h 

If the values of l\ and I 2 are both unknown, let l\(n) be the numerical approximation 
of the solution of the following equation : 


( 3 ) 


In - 2/1 1 1 

yi~h(n) 2 2 Xp 


2M„ 


In - k{n)_ 


= 0. 
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It defines an approximation of l\ because Equation ([3]) is a discretization of Equation 
([TJ. Moreover, by Corollary 2.3 : 

2S 

h(n ) := 2'y n -li(ri), m n := --— ^ " and p n := log 


In - h(n) 


yi~h{n) 2 [ 7 „ — h(ji)\ 


define some converging approximations of I 2 , m and p respectively. 


3. A MODEL OF THE DOSE-EFFECT RELATION 

The model introduced in this section is tailor-made to study the relation between 
a dose d > 0 and a quantitive measure E(d) of an effect of the administered drug. 
Precisely, E(d) can be a measure of the therapeutic response or a measure of a 
side-effect. 


3.1. The skew normal distribution. On the skew normal distribution, see Az- 
zalini [2] and Chen et al. [5]. 

Definition 3.1. The skew normal distribution of parameters f £ M (location), 
to > 0 (scale) and a £ ffi. (shape) is the probability measure II on (M, S(K)) defined 
by : 


n(dx)= 


dx 


CN 

1 + erf 

'a(x-C)!' 

V 2 w 2 J 

x>\J 2 



t 

The skew normal distribution of parameters £, oj and a is denoted by Af(f, w 2 , a). 

Consider £, a £ R. and cu > 0. The parameters of the skew normal distribution 
Af(£,,uj 2 , a) are related to its mean p, its standard deviation cr and to its skewness 
coefficient 7 (Pearson) as follow : 


S 2 v 2 A c x 

a = , , to = - - „. and F = p — uio 

v / T^ 2 1 — 2 A 2 /tt 


with 


\S\ ■= 


l7| 1 / 3 %/V2 


\/|7| 2/3 + [(4 — 7t)/2] 2 / 3 ' 

Let (Xi,..., X n ) be a n-sample (n £ N*) such that X\ A /"(£, uj 2 , a). Consider 


^ 1 


X n := - V X t , S 2 n := - V(Xi - X n ) 2 and % := - V 

n L ' n z ' n z ' 


Xi-x r : 3 

Sn 


So, 


OLn • — 


1 - 6) 


D “I ~ 


Sl 


with 


\S n \ := 


1 - 25 2 /tt 

| 7 n| 1 / 3 \AV 2 


and : — X n ca n S n ^j ^ 


V'l7n| 2 / 3 + [(4-TT)/2] 2 /3 

are some converging estimators of a, w 2 and £ respectively. 


3.2. The skew normal model of the D-E relation. The model is a family of 
random variables E := {E (d),d £ R + } satisfying the following assumption : 

Assumption 3.2. 

(1) For every d £ R + , E(d) A f(f(d), uj 2 (d), a{d)) where £,o £ C°(M_|_,R) 

and uj £ C°(R+, ]0, 00 [). 

(2) p : d £ M+ 1 — > E[E(d)] is a logistic function. 
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( 3 ) 


There exists do £ R + such that a : d £ R + i —> y/var[E(d)] is decreasing on 
[do, oo [, and 

lim a (d) = 0 . 

d—> oo 


Assumption 3.2 (3) means that the variability of the measured effect between 
the patients tends to vanish for large doses of the drug. 


Assume that E(d) is a measure of the therapeutic response of the administered 
drug at the dose d £ R+. It is usual to assume that E (d) is a Gaussian random 
variable. The skewness of the empirical distributions are taken into account in the 
model E of the therapeutic response in order to refine the choice of an optimal 
dose. Indeed, the skewness of the distribution of E(d) indicates if the therapeutic 
responses of the major part of the patients are over (positive skewness coefficient) 
or under (negative skewness coefficient) the mean therapeutic response. Therefore, 
if the mean of E(d) is high enough, its skewness coefficient is positive and its stan¬ 
dard deviation is small enough, then the dose d is admissible. 


With the notations of Assumption |3.2| the last part of the section deals with the 
estimation of the functions /i and er, and then of the functions £, u> and a. So, 
for each admissible dose d > 0, Assumption 3.2 (1) allows to simulate the measure 
of the therapeutic response E(d) multiple times for a better evaluation at the dose d. 


Consider the finite set D C R+ of the administered doses and, for every d £ D, 
n £ N* independent observations ei(d),..., e n (d) of the random variable E(d). For 
each dose d £ D, m(d), <r(d) and 7 (d) can be estimated by 


respectively. 


m„(d) 


£ n («0 

7 n(d) 


n 

^c(d), 

i= 1 

1 n 

- ^[e,(d) - fhn(d)] 2 and 
U i= 1 

ej(d) - m n (d) l 3 
CT n (d) 



Since m is a logistic function by Assumption |3.2| its parameters can be approx¬ 
imated by the method stated at Section 2. 


Assume there exists do £ D such that 

d £ D 1 —> a n (d) 

is decreasing on D n [d 0 , 00 [. 

(1) If the function d £ D 1 —> <7„(d) is constant on D (~l [0,do[, then a could 
be a logistic function. Its parameters can be approximated by the method 
stated at Section 2. 

(2) If the function d £ D 1 —> a n (d) is increasing on D D [0, do[, then a could be 
a Gaussian-type function : 

a (d) := l + exp(—md 2 + pd + q) ; Vd € M+ 
with m > 0 and l,p,q £ R. Since 

lim u(d) = 0 
d—> 00 
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by Assumption 3.2 (4), then 1 = 0. So, the parameters to, p and q can be 
approximated by polynomial regression on {(d,y n (d)); d G D}, where 

y n {d) := \og[d n (d)] ; Vd € D. 

This method is classic (see P. Armitage et al. [T|). 


4. An APPLICATION TO AN ANTIDEPRESSANT MEDICATION 

This section deals with a model of the escape time (ET) of rats under an exper¬ 
imental antidepressant (AD) medication. For confidentiality reasons, the name of 
the drug is not specified. 


The escape time (in seconds) is interpreted as a therapeutic response to the AD. 
The AD is administered at four different doses : 0, 0.75, 1.5 and 3 mg.kg -1 . The 
escape time is evaluated on 32 rats ; 8 rats for each posology. 


For each dose, the mean, the standard deviation and the skewness coefficient of 
the escape time have been computed : 


Statistics Doses 

0 

0.75 

1.5 

3 

Mean 

33.3875 

44.1625 

51.5 

78.225 

Standard deviation 

26.9715 

30.8113 

44.6582 

31.9657 

Skewness 

-0.0276 

-0.1381 

1.2827 

0.3504 


4.1. The model of the dose-escape time relation. According to Assumption 
|3.2[ the mean escape time is modeled by a logistic function of the administered 
dose. By using the approximation procedure stated at Section 2 : 

TOET(d) = 21.8153 + 0 0116 + e — 0 . 8278 ^— 2.5929 ? W € R +‘ 

The standard deviation and the skewness coefficient are modeled by Gaussian type 
functions of the administered dose. By using polynomial regressions as suggested 
at Section 3 : 

cr ET (d) = exp(—0.1502 • d 2 + 0.5289 • d + 3.2459) 

and 

7 ET (d) = 0.2381 + exp(—0.6503 • d 2 + 2.2935 • d - 1.5578) 
for every d G R+. 

The functions TOet, <tet an d 7et are Pitted on the interval of doses [0,4] : 



Figure 1. Mean ET 
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Figure 2. Std. dev. ET 


FIGURE 3. Skewness ET 


4.2. The optimal dose. The optimal dose should maximise the mean and the 
skewness coefficient, and minimize the standard deviation of the therapeutic re¬ 
sponse. 

Unfortunately, there is no dose such that these conditions are satisfied together 
by the model of the dose-escape time relation. 

Consider the dose d opt := 3 which maximizes the mean on [0,3], and such that 
^ET(^opt) = 32.4903 and 7 E T(d op t) = 0.3504. 


Since the AD tested in the clinical trial has a priori no significant side-effect for 
the doses less or equal than 3 mg.kg -1 , cZ opt could be the optimal dose. Indeed, it 
maximizes the mean of the escape time, aET(rfopt) c l° ser to 

min <tet(gO = 26.9715 

de [o,3] 

than to 

max <tet(gO = 44.6582, 
de[0,3] 

and 7ET(rf op t) > 0. 


Appendix A. The proofs of Section 2 
Proof of Proposition \ 2. 2\ For every x £ R, 

me mx+p 


f'{x) = - 


and 


/"(*) = ] 


[(/ 2 - h)- 1 + e mx +P} 2 


1 


a mx+p 


[(k ~ h)- 1 + e mx +P} 3 \l 2 - h 
(1) Assume that m < 0. So, 

lim f(x) = h + 0 = l\ 


J 


and 


lim f{x) =h + - - 

x—¥oo (y 2 — Ill 


Assume that m > 0. So, 

lim f(x) = l\ + 


x—>—oo 


(h — h) 1 


= h- 


= I2 














A SKEW NORMAL MODEL OF THE DOSE-EFFECT RELATION IN PHARMACOLOGY 9 


and 

lim f(x) = li + 0 = l\. 

x—too 

(2) If to < 0 (resp. to > 0), f'{x) > 0 (resp. /'( x) < 0) for every x £ R. So, if 
to < 0 (resp. to > 0), then / is increasing (resp. decreasing) on R. 

(3) f"{x) = 0 if and only if x = 9. Moreover, if to < 0 (resp. to > 0), f"{x ) > 0 
if and only if x s] — oo, 9[ (resp. x g ]9, oo [). So, the graph of the function 
/ has a unique inflection point, at 9. 

That achieves the proof. □ 


Proof of Corollary \2.3\ By Definition |2.1 

(4) /(0) = l 


By Proposition 2.2 

(5) 

( 6 ) 

(7) 


1 (h-h^+eP- 

\og(l 2 -h)+p 


f(9) = and 

/'(<?) = m ^ 2 ~ Z i) 


By Equation Q, l 2 = 2/(0) — l-\. By Equation ([7| : 

0 )*=* m = ~m 9 )-h] 


m = — - 


2 

2 f(9) 


m-h ' 

By Equation Q, since l 2 > /(0) by the definition of / : 

1 


@ 


/( 0 ) = h + 

p = log 


1/2 [f(9)-h]-i+eP 

1 1 


By Equation 0 : 


m-h 2[f(9)-h]\' 


0 


1 


9 = -log 

TO 


o /w-zi ; 

. /(o) - h 


log 


o /( 0 w(°),; 

/(0 )-(i + 


2 9f(9) 

m - h 


= 0 . 


That achieves the proof. □ 


Proof of Proposition \2.J\ For every x € R, 


f(x) = - 


mx+p 


[(l 2 - h )- 1 +e mx +P ] 2 


= —To 


h + 


1 


(l 2 - h)- 1 + e mx+p 


-h 


1 - 


1 


l 2 — h 
= —m[f(x) - h] 


h + 


1 


1 - 


(h - ^ i )' 1 + e mx +P 
1 


-h 


l 2 — h 


[f( x ) ~ h 


So, the logistic function / is the solution of Equation 0. □ 
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